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Numerical simulations of a dissipative hard sphere gas reveal a dependence of the cooling 
rate on correlation of the particle velocities due to inelastic collisions. We propose 
a coefficient which characterizes the velocity correlations in the two-particle velocity 
distribution function and express the temperature decay rate in terms of this coefficient. 
The analytical results are compared with numerics. 

Keywords: kinetic gas theory, many-particle systems, dissipative gases, granular gases, 
velocity correlations 

1. Introduction 

The Boltzmann equation for the velocity distribution / (v, t) of particles of a gas has 
been derived with the use of the molecular chaos assumption, which implies that 
all correlations between the velocities of particles may be neglected 1 . This equation 
was also a starting point to build up a kinetic theory of dissipative gases where the 
particles collide inelastically (for an overview see, e.g. 2 ). Whereas the molecular 
chaos assumption is justified for dilute molecular gases, there are two reasons for 
its violation for the case of moderately dense dissipative gases: First the positions 
of the particles are not independent due to excluded volume effects which gives rise 
to the Enskog factor. 3 A second source of correlations is the dissipative nature of 
particle interaction. Since dissipation concerns only the normal component of the 
relative velocity of colliding particles, the angle between their velocity vectors after 
the collision is always smaller than the angle just before the collision. The latter 
effect leads to aligning of the trajectories of adjacent particles which contradicts 
the molecular chaos assumption. Obviously this effect increases with increasing 
inelasticity, expressed by the coefficient of restitution e which relates the normal 
component of the relative particle velocity after and before a collision. This expec- 
tation has been confirmed by studies of randomly driven Granular Gases 4 , where it 
was shown that the molecular chaos assumption becomes less accurate with increas- 
ing degree of particle inelasticity. Instead, the velocities of the colliding pair are 
noticeably correlated; this effect depends on the restitution coefficient and on the 
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gas density. 5 On the macroscopic (hydrodynamic) scale these velocity correlations 
cause vortex structures in the velocity field of a spatially homogeneous dissipativc 
gas. 6 

In spite that obviously the microscopic correlations of the velocities of adjacent 
particles play an important role for the kinetics of dissipative gases, the particu- 
lar form of the multi-particle velocity distribution function, which describes these 
correlation is not known. It is the aim of this paper to present the first (to our 
knowledge) attempt to construct the two-particle correlation function of a simple 
dissipative gas. 

2. Boltzmann-Enskog equation for hard spheres 

In kinetic gas theory the evolution of the velocity distribution function f (v,t) is 
governed by the Boltzmann equation 1 

= /(/,/) (1) 

with the collision integral 

/(/, /) = a 2 J dv 2 J deO (-v 12 ■ e) \v 12 ■ e\ [ X f K, t) f (v' 2 \ t) - f (v u t) f (v 2 ,t)} . 

(2) 

The velocities v\ and v 2 are the pre-collision velocities of the direct collision, and 
the unit vector e points from the center of particle 2 with velocity v 2 to the center 
of particle 1 with velocity v\. Here v\ 2 = v\ — v 2 , a is the particle diameter and 
Q(x) is the Heaviside step-function. The velocities v", v 2 correspond to pre-collision 
velocities of the inverse collision, which ends up after the collision with the velocities 
vi and v 2 . 

For a gas of inelastically colliding particles the function 

= \g"\ V{v'{,^) 
X ~ \g\ V(v,,v 2 ) W 

accounts for the Jacobian of the transition from the pre-collision velocities v'{, v^' 
of the inverse collision to those of the direct collision v\ , v 2 and for the ratio of the 
lengths of the collision cylinders of the inverse and the direct collision. These are 
proportional to the normal components of the relative velocities, i.e. g = (v~i 2 ■ e) e 
for the direct collision, with the analogous expression for the inverse collision. The 
ratio 

e=\g\/\g"\ (4) 

defines the coefficient of restitution e which describes the inelasticity of the particles. 
Using this coefficient the velocities of two colliding particles before and after a 
collision can be related. In particular, these relations read for the inverse collision 

v*l = vi - {v 12 ■ e) e (5) 

4' = #2 + {vi2 ■ e) e. (6) 
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The coefficient of restitution is the central quantity in the theory of granular gases. 
Experiments 7 as well as theory 8 ' 9 * 10 show clearly that the coefficient of restitution 
is not a constant but depends on the impact velocity. Nevertheless, for the sake of 
simplicity of calculation it is frequently assumed that this coefficient is a material 
constant. With the assumption e = const, from (3) we find \ = e ~ 2 - 



3. Correlations in dissipative gases 

The most important precondition for the derivation of (1,2) is the assumption of 
molecular chaos (Stofizahlansatz) which expresses the two-particle distribution func- 
tion as a product of the one-particle distribution functions 

f2(vi,V 2 ,t) = f(v 1 ,t)f(v 2 ,t) (7) 

which is identical with the assumption of vanishing correlations between successive 
collisions of each particle. For gases of particles of finite sizes this assumption is 
violated due to finite volume effects (i.e. two particles can not occupy the same 
volume). This gives rise to the Enskog factor 11 (72(c) 

/ 2 (tfi,tfc)« 92(a) f (V!,t) f (v 2 ,t) (8) 

which accounts for the increasing collision frequency in such gases. This quantity 
equals the contact value of the pair correlation function 32(0"). For systems of hard 
discs it reads: 12 



where r\ = imo 2 /4 is the packing fraction, n = N/V is the number density with iV 
being the number of particles in the system and V is the two-dimensional volume 
(area). The Enskog factor obviously accounts only for static spatial correlations. 
Any dynamic correlations between the velocities of particles are still ignored. 

For the case of dissipative gases there is an additional source of correlations. 
These arise since the restitution coefficient affects only the normal component of 
the relative velocity of a colliding pair, leaving the tangential relative velocity un- 
changed. If two inelastic particles collide with a (pre-collision) angle a" the angle 
a after the collision will be smaller, a < a" (see Fig. 1). 

As a consequence subsequent collisions lead to more and more aligned traces 
of neighboring particles. It has been shown 6,13 that these correlations give rise 
to convection cells in an otherwise homogeneous granular gas, i.e., to macroscopic 
correlations in the velocity field even in absence of spatial inhomogeneities. 

Since the particles are aligned in macroscopic vortices, such velocity correlations 
lead to a reduced collision frequency as compared with the collision frequency of 
a gas of elastic particles of the same temperature. The reduction of the collision 
frequency leads in its turn to a reduction of the cooling rate, i.e., to a retarded decay 
of temperature even in the stage of homogeneous cooling. Therefore, to estimate 
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Fig. 1. Elastic (left) and inelastic (right) spheres before (black), during (dark grey) and after a 
collision (light grey). For elastic spheres the incoming angle equals the outgoing angle (a = a") 
whereas for inelastic collisions a < a" . 



the importance of velocity correlations we analyze the cooling rate of a dissipativc 
gas. 

4. Temperature decay of dissipative gases 

Assume a homogeneous dissipative gas of particles with velocities which obey a 
distribution function / (v). For e = const, the Boltzmann equation decouples into 
a set of simpler equations 14 , one for the time dependent temperature 



m f 
= 3W 



T(t)=— / dvv 2 f(v,t) (10) 



and another one for the time independent (normalized) distribution function. With 
the assumption that the velocities obey a Maxwell distribution it has been shown 
by scaling arguments that temperature decays due to dissipative collisions as 

with To being the initial temperature and To being a scaling factor which will be 
discussed in detail below. Equation (11) is known as Haff's law 15 . Later deviations 
from the Maxwell distribution have been found 16,17 , however, the temperature decay 
law is still correct. 

Let us compare the analytical results for the temperature evolution with the 
temperature decay as obtained by an event driven Molecular Dynamics simulation. 
Haff's law was derived with the molecular chaos assumption, therefore, in order to 
obtain results which are comparable with theory we apply a special numerical proce- 
dure which destroys the collision-induced velocity correlation of colliding particles. 
This procedure is illustrated in Fig. 2. We shuffle the velocities of the colliding par- 
ticles with the velocities of randomly selected particles. It is easy to understand that 
the shuffling procedure is unique since the condition for the after-collision normal 
relative velocity to be positive must be fulfilled. 



VIOLATION OF MOLECULAR CHAOS IN DISSIPATIVE GASES 























Fig. 2. The velocity shuffling procedure. Left: the black particles at the instant of the collision 
with pre-collision velocities. Middle: the colliding particles are assigned new velocities due to the 
collision rule Eq. (5), two randomly selected particles arc drawn dark grey. Right: the velocities 
of the colliding particles are shuffled with the velocities of the randomly selected particles. The 
shuffling procedure is unique due to the condition for the after-collision normal relative velocity 
to be positive. 




Fig. 3. Evolution of the temperature. Dashed line: event driven MD simulation (e = 0.93, 
n = 0.05), full line: Haff's law, Eq. (11) together with the randomized event driven simulation 
(lines fall together). As soon as clusters appear (hit fa 6) there is a large deviation of the numerical 
results from theory. The inset shows a magnification of the range 2.5 < lni < 3. In the entire 
range of time we observe a small but systematic deviation between the curves. 
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The solid line in Fig. 3 shows the evolution of temperature for a system of 
N = 5, 000 inelastically colliding spheres in two dimensions with periodic boundary 
conditions where the described shuffling procedure has been applied (randomized 
simulation). The reduced density was n = Na 2 /V = 0.05 (with a = 1), where V is 
the area of the quadratic box and the coefficient of restitution is e = 0.93. In the 
entire range of t the curve is identical up to the line width with Haff's law (11) where 
the relaxation time To is a fit parameter. Repeating the same simulation without the 
shuffling procedure we obtain a temperature decay as drawn by the dashed line in 
Fig. 3. The curves are close to each other until \nt s=s 6. At this time instant for the 
unshufflcd system we observe formation of spatial inhomogeneities (clusters) which 
leads to a reduced collision frequency. The cluster formation in dissipative gases 
due to a pressure instability is a well understood phenomenon which was described 
first in 18 . 

Obviously, our randomization procedure suppresses the formation of density 
inhomogeneities which can be understood from mode-coupling analysis. 18 We know 
that the space inhomogeneities arise in two steps: first vortex patterns appear 
due to interparticle velocity correlations induced by the collisions, then clusters 
develop owing to interaction between the density modes and the flow- velocity modes. 
Thus, the collision-induced correlations play a clue role in the formation of spatial 
inhomogeneities. Therefore, since clustering or vertex formation is not possible 
in the randomized simulation due to suppression of the collision-induced velocity 
correlations we find almost perfect agreement with Haff's law. 

Figure 3 shows that in the region hit < 6 Haff's law describes the temperature 
decay of the dissipative gas with and without randomization up to a good degree 
of accuracy. Having a closer look to Fig. 3 (inset) one recognizes, however, a 
small but systematic deviation between these two temperature curves. We have 
checked spatial homogeneity of the particle distribution for the both gases and did 
not observe any difference in the time span where Haff's law is applicable. Thus we 
conclude that this difference may be attributed only to velocity correlations of the 
particles of the gas. 

It is possible to fit the value of the relaxation time To in Haff's law to reproduce 
the simulation data (dashed line) up to the same accuracy as done before for the 
randomized gas. The numerical value of To which depends on the coefficient of 
restitution and on the density of the gas may, therefore, serve to characterize the 
velocity correlations of the particles. We will return to this issue in Sec. 6. 

5. Beyond the concept of molecular chaos 

To develop a theory beyond the concept of molecular chaos we introduce the reduced 
velocities 

c= — with vo = 

v 

where v is the thermal velocity and the reduced velocity distribution function / (c) 
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which is defined by 



n 



f(^) = - m f(c,t) . (13) 
v o\ l ) 

In the equivalent way we define the reduced two-particle distribution function 
h (ci, c 2 , ri2, t) with fi2 = fi - f 2 . 

Correlations of the particle velocities are taken into account by assuming the 
following form of the two-particle velocity distribution function: 

f2(ci,c 2 ,r 12 ,t) = f(c 1 )f(c 2 ) 92 (ri 2 ) +x( r i2)D (r 12 ^ : cic 2 (14) 

The first term in the brackets in the right-hand side of Eq. (14) accounts only for 
static correlations, while the second term accounts for dynamic correlation of a pair 
of particles with the reduced velocities ci and c 2 . The traceless dyad of the unit 
vector r\2 — r\2jr\2 is defined as 



d (n 2 ) 



^r 12 ri2-U (15) 



with U being the unit tensor and we use the conventional notation of vector algebra, 
i.e., aa denotes the dyad and 

D (?i2^ ■■ c*ic 2 = 3 (ci ■ Fi 2 ) (c 2 ■ f 12 ^ - (ci • c 2 ) . (16) 

This is one of the simplest forms which may be built up from three vectors and which 
features the requested symmetry with respect to vectors ci and c 2 . The function 
X (T12) expresses the dependence of the velocity correlations on the particle distance 
which, obviously, vanishes for infinite distance, i.e., x( r i2) — * for T\ 2 — > 00. 
At the instant of a collision, i.e., when r\ 2 — (or ri 2 = e) Eq. (14) reads 

k (ci, c 2 , ae) = / (ci) /(c 2 ) g 2 (a) {1 + b [3 (5i ■ e) (c 2 • e) - (ci • c 2 )]} , (17) 

where 6 = x (f) /ff2 (o - )- For the one-particle distribution function we use its repre- 
sentation in terms of the Sonine polynomial expansion with the first non-vanishing 
Soninc coefficient a 2 : 

/(ci)=</»(c)[l + a 2 S 2 (c 2 )] . ( 18 ) 
where for the two-dimensional case the second Sonine polynomial reads, 

S 2 (x) = ^x 2 -2x + l. (19) 

With the distribution function given by Eqs. (17,18,19) all necessary quantities 
which describe the kinetic properties of the gas can be obtained, e.g., the second 
moment of the collision integral 

M2 = ~ f dci J dc 2 J deO (-c\ 2 • e) |ci 2 • e\ f 2 (c 1 ,c 2 ,ae) A (c\ + c 2 ) , (20) 
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where A (c\ + c 2 ) denotes the change of the quantity c\ + c 2 at the dissipative 
collision. 

The importance of the second moment yu 2 follows from the fact that it charac- 
terizes the dissipation rate due to inelastic collisions. Indeed, the factor 9 (— c*i 2 • e) 
in the integrand guarantees that particles with velocities c\ and c 2 do collide. The 
factor |ci2 • e| characterizes the length of the collision cylinder, which multiplied by 
the factor / 2 (<?i , 62 , ere) characterizes the number of collisions between particles with 
velocities Ci, <?2 and the intercenter vector ere, occurring per unit time. The last 
factor in the integrand A (c\ + c 2 ) characterizes the energy loss in such collisions. 
Integration over all possible velocities and intercenter vectors yields the quantity 
which characterizes the energy loss per unit time, i.e., the decay rate of the gas 
temperature. 

Straightforward calculations yield 



1 3 63 1 

1+ 16 a2 -40 6 



(21) 



This expression shows that ^2 depends on both coefficients 02 and & in a similar 
way. 

6. Influence of the velocity correlations on the temperature evolution 

Originally, Haff's law (11) was derived by scaling arguments with the assumption 
of the Maxwell distribution for the particle velocities. 15 From this derivation it is 
not possible to obtain the relaxation time To in terms of the microscopic particle 
properties. Starting with the Boltzmann-Enskog equation, however, the tempera- 
ture evolution 17 , i.e. Haff's law, and expression for the relaxation time To in terms 
of the second moment of the collision integral can be found: 



3 ' V m 
With expression (21) for [i 2 we obtain 

1 

1-e 

with 



- 1 = i 32 (a)a 2 ^/^° M2 . (22) 



to = :r— 2 A{e) , (23) 



A- 1 (e)=uj a 



3 63, 
1 + T6 a2(£) -40 5(e) 



The second Sonine coefficient for a two-dimensional gas of inelastic discs as a func- 
tion of the coefficient of restitution has been obtained with the molecular chaos 
assumption (i.e. ignoring particle correlations) 16 : 

16(1 -6) (1-26*) 

a2 ^ 57- 256 + 3Ck2(l- e )- (25) 
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Using Eqs. (23,24,25) the (macroscopic) relaxation time To is expressed in terms of 
microscopic quantities which characterize the particle collisions. 

Analyzing Eqs. (24,25) it turns out that the dependence of A(e) on e is very 
weak: In the interval e e (0.8, 1.0) it varies within less than 0.6%. From Eq. (24) 
follows that the velocity correlations, expressed by b preserve the functional form 
of the temperature decay, but change the value of the relaxation time. Moreover, it 
also follows that the more pronounced the velocity correlations (i.e. the larger the 
coefficient b) the more the relaxation time To differs from the relaxation time of a 
gas where the velocity correlations are suppressed. Taking into account that <Z2 is 
small and that b is also not large, we can write for the factor A(e): 



3 , . 63 , . . 
16 fl2(£)+ 40 fe(£) 



= Aand + 1^(7 ( 2 6) 



where A ranc j is the factor A that corresponds to a gas without velocity correlations, 
i.e., a randomized gas. 

This theoretical prediction can be checked by numerical simulations: again we 
perform event driven Molecular Dynamics simulations, measure the temperature as 
a function of time and determine the value of t as a fit parameter due to Eq. (11). 
According to Eq. (23) the function A(e) can be obtained by drawing To (l — e 2 ) 
over e. The results are shown Fig. 4 for two values of the density n. 
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Fig. 4. The function A(e) = To (l — e 2 ) over the coefficient of restitution e for a simulation of 
5,000 inelastic hard disks of density n = 0.05 (left) and n = 0.1 (right). The lines have been fit 
to the numerical data. For the randomized simulation (squares) the curve does not significantly 
depend on e whereas in agreement with the theoretical prediction Eq. (26) the original gas shows 
noticeable dependence on e due to velocity correlations. 



The squares corresponds to the randomized gas and do not exhibit any notice- 
able dependence on the coefficient of restitution e. As it has been already mentioned 
this is a consequence of the very weak dependence of A on e via the second Soninc 
coefficient 02 . The triangles show the same data for simulations where the random- 
izing procedure has not been applied. The lines show linear least square fits to the 
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numerical data. Figure 4 shows that the dependence of A on e via the correlation 
coefficient b can be approximated by a linear function of e. 

7. Conclusions 

We study microscopic particle velocity correlations in dissipative gases. Taking 
into account velocity correlations preserves the functional form of the temperature 
decay, which still obeys Haff's law but changes the temperature relaxation time. 
Going beyond the assumption of molecular chaos we suggest the functional form 
of the two-particle distribution function which accounts for velocity correlations 
and calculate the temperature relaxation time for this distribution function. We 
performed Molecular Dynamics simulations of dissipative gases and suggest a ran- 
domizing procedure which suppresses the velocity correlations but leaves all other 
properties of the gas dynamics unchanged. From the comparison of the results of 
the simulation of the randomized gas with traditional simulations we obtain the 
dependence of the correlation coefficient on the coefficient of restitution. 
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